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D , Group field theory is a generalization of matrix models, with triangulated pseudomanifolds as Feyn- 

' man diagrams and state sum invariants as Feynman amplitudes. In this paper, we consider Boula- 

tov's three-dimensional model and its Freidel-Louapre positive regularization (hereafter the BFL 
model) with a 'ultraviolet' cutoff, and study rigorously their scaling behavior in the large cutoff 
limit. We prove an optimal bound on large order Feynman amplitudes, which shows that the BFL 
model is perturbatively more divergent than the former. We then upgrade this result to the con- 
structive level, using, in a self-contained way, the modern tools of constructive field theory: we 
construct the Borel sum of the BFL perturbative series via a convergent 'cactus' expansion, and 
^ 1' establish the 'ultraviolet' scaling of its Borel radius. Our method shows how the 'sum over trian- 

U I gulations' in quantum gravity can be tamed rigorously, and paves the way for the renormalization 

. program in group field theory. 
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I. INTRODUCTION 



(N 
> 

■ A. Renormalization in quantum gravity 

. 

• ' Since the quantum theory of metric perturbations about Minkowski space-time was proved non-renormalizable in 
four dimensions research in quantum gravity has followed divergent avenues, guided by different reactions to 
this frustrating no-go result Q]- Supergravity and string theory have relied on extended symmetries and Kaluza- 
Klein dimensions to tame the ultraviolet (UV) behavior of gravity, even perhaps subsuming it, together with all 
• • r other particles and interactions, under the unifying concept of strings 0, Non-commutative geometry has explored 
. ^ ' generalizations of space-time in which UV divergences might disappear, and in which the standard model becomes more 
^ , natural [1, Q. The putative existence of a non-Gaussian UV fixed point of quantum general relativity (GR), coined 
5^ ' the "asymptotic safety scenario" by Weinberg, has been explored in computer simulations, with intriguing results 
5^ , Loop quantum gravity (LQG) has followed the historical path, opened by Dirac, of the canonical quantization of the 
full gravitational field, discarding the issue of renormalizability together with the Minkowski back grou nd [1, [l^] ■ In 
other background-independent approaches, such as causal dynamical triangulations [ll| or causets the framework 
of quantum field theory (QFT) is abandoned altogether, and space-time is envisaged as a combinatorial build-up, in 
the spirit of condensed matter physics. 

Even if none of these approaches can claim to have reached maturity, it is already clear that deep and fruitful 
advances have occurred along the way. Let us mention but a few of them. On the mathematical side, the development 
of topological quantum field theory has led to major discoveries in the theory of manifolds and knot invariants, initiated 
by il3|. On the physical side, the old concept of 'space-time discreteness' has become much more precise, to the point 
where the spectra of area and volume operators are explicitly computable [13, HI, [ll]; the non-commutative field 
theories ^IT'I inspired from string theory and noncommutative geometry have 'exorcised' the Landau ghost fis', T^; 
the discovery of the AdS/CFT correspondence has provided a fresh new look on strongly-coupled gauge theory |20j; 
convincing scenarii for the quantum resolution of space-time singularities have been put forward [2l|. And so forth. 

In this impressionistic picture of perspectives on quantum gravity, however, a pattern strikes us: the less they 
rely on background geometric or topological structures, the more remote from standard field-theoretic ideas they are. 
In particular, the concept of the renormalization group (RG), that is of the hierarchical integration of fluctuations 
about a computable solution, seems to have somehow faded away.^ This Wilsonian paradigm, however, has proved 



Obviously, this is not the case in the asymptotic safety scenario. The recourse to uncontrolled truncations of the effective action, 
however, makes this program difficult to justify from a mathematical physics perspective. 
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fundamental in many branches of theoretical physics, and is arguably the only encompassing framework to address a 
system with many interacting degrees of freedom. 

Can we find a field-theoretic formulation of quantum gravity, free of a priori geometric or topological background 
structure, in which a RG flow could be rigorously defined? In this paper, we take some first steps towards this 
goal. Our starting point, group field theory (GFT) [13, li^s at the crossroads of matrix models (and dynamical 
triangulations), and LQG, in its spin-foam formulation. 

B. Group field theory 

Quantum gravity should implement a sum over all geometries of space-time with certain weights. From this 
perspective, general relativity is only of limited guidance: although the Einstein-Hilbert action suggests how to weigh 
different metrics on a fixed differential space-time manifold, it does not tell us how to weigh the sum over different 
topologies of space-time (or over differential structures in dimension 4 or above) . This problem is loosely called the 
third quantization of gravity: quantum gravity may not only be about quantizing metric fluctuations over a classical 
background topology; perhaps it should also quantize the space-time topology itself. 

The discovery that matrix, not scalar models, have a topological content 1241 opened a window on how to solve this 
third quantization problem with field-theoretic methods in two dimensions [25| . In matrix models, Feynman graphs 
are ribbon graphs, dual to triangulations of surfaces. As the size of the matrix increases, the scaling of observables 
favors the particular class of planar graphs which correspond to triangulations of the sphere, higher-genus surfaces 
being exponentially suppressed. 

The natural candidate generalizations of matrix models in higher dimensions are tensor models, and we propose to 
nickname QGd the would-be d-dimensional model implementing the third quantization of gravity. QGd should follow 
the rationale of matrix models, namely that space is fruitfully described by gluing together more elementary "cells" 
through attaching maps. In the combinatorially simplest case, these cells are just d-simplices. Since a d-simplex has 
d+1 facets on its boundary, the backbone of our would-be QG^ should be some abstract (/)'^+^ interaction on rank-d 
tensor fields (j). Now, the GFT's of Boulatov and Ooguri [2^, [23| provide ansatze to generate precisely such tensor 
fields: (/> is a scalar field over d copies of some compact group, and harmonic analysis yields rank-d tensors, which can 
be naturally coupled in the combinatorial pattern of d-simplices. Of course, the exact forms of the propagator and 
interaction in this setting remain to be found. One should also include a finite list of relevant and marginal terms to 
ensure renormalizability and probably constructive stability (sec below). Ideally, these new terms would also imply 
(or at least allow and favor) the standard model matter content of our universe. 

Surprisingly enough, LQG, which takes a radically different perspective on quantum gravity, notably by keeping a 
classical background topology, has converged with the GFT framework. In LQG, the kinematical states of the quantum 
gravitational field are well understood in terms of spin networks, and the main open issue is the implementation of 
dynamics, through the Hamiltonian constraint. Spin- foam models provide heuristic definitions of such dynamical 
transition amplitudes, obtained via a discretization of general relativity on a triangulation of space-time (28l . [29| . 
There are now interesting such spin-foam models [13, [U, [12, [13 j hereafter called EPR-FK models, which in four 
dimension reproduce Regge gravity in a certain semiclassical limit (s^ . [35| . There are also some glimpses that they 
might be just renormalizable [36| . These spin- foam models, however, capture only a finite subset of the gravitational 
degrees of freedom, and the question arises of the existence of a 'continuum limit'. As the spin-foam amplitude can 
always be interpreted as a Feynman amplitude of a suitable GFT |37i] . this question boils down to the problem of 
renormalization in GFT. 

This is what makes GFT's exciting: their Feynman amplitudes are spin-foams, which implement a tentative 
background-independent second quantization of gravity, and they are scalar field theories of the right tensorial type 
with the right 0^^+^ interaction for third quantization! They thus seem ideally suited for performing both second and 
third quantization of gravity in a single bold move. This coincidence also leads to an even bolder hope, namely that 
perhaps adding third quantization is just the missing ingredient to the second quantization problem as well. 

C. Multiscale renormalization 

To begin the program of renormalization of GFT, the first step is the definition of scales. However, because third 
quantization blurs space-time, it fools our deeply ingrained relation of scale to distances (centimeter scale, millimeter 
scale...). In particular our familiar intuition that ultraviolet is about short distance and infrared about long distance 
effects may no longer hold. Fortunately non-commutative field theory, where ultraviolet-infrared mixing can occur, 
as in the Grosse-Wulkenhaar model [38], has familiariazed us with such conceptual issues. It provides a hint as to the 
extended definition of scales and renormalization group (RG) which are presumably required for QGd |39l|. 
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Roughly speaking, to set up a RG analysis in QFT, one must draw a line between propagator and interaction. 
There should be no scales in the interaction, because vertices should always be considered, by definition, the correct 
"local" objects (even when, like in the x-space representation of the GW model, they naively don't look local at all). 
Hence scales should be entirely contained in the propagator, more precisely in its spectrum, which does not depend on 
the basis used. In standard Euclidean QFT, a propagator is a positive covariance C which is the inverse of a positive 
Hermitian operator iJ, and we are led to define a physical scale as a range of eigenvalues for H. A RG trajectory 
should then integrate over 'ultraviolet' scales (large eigenvalues) towards 'infrared' scales (smaller eigenvalues). The 
best way to cut the typical geometrically growing sequence of RG slices is through the parametric representation of 
the propagator as 

C = H-^= / e-^-f^da = VC,,Co = / e'^^da, C, = / e'^'^da for i>l, (I.l) 

Jo Ji Jm-^' 

where M is a fixed number greater than 1^. Ci should be considered the definition of a (discrete) RG scale. The 
RG trajectory then performs functional integration over scales, one at a time. A renormalizable trajectory is one in 
which the relevant part of the RG flow changes only a finite number of terms. In beautiful cases like QCD, such a 
trajectory joins a perturbatively computable UV fixed point (asymptotic freedom) to a non-trivial phase transition 
(quark confinement). We hope for such a scenario, namely that a consistent QG4 flow could join some UV fixed point 
(presumably topological or at least with enhanced symmetry) to a non-trivial phase transition at which our ordinary 
effective geometry would emerge.^ 

Within this program, we hope to clarify the meaning of renormalizability for third quantized theories, and in 
particular to answer whether or not the recent EPR-FK models define a renormalizable quantum theory of gravity, 
to understand how they should be physically interpreted, and eventually to find out whether and how they should be 
modified for full mathematical consistency up to the constructive level. 

This is an ambitious program and it is reasonable to gain better expertise first on QG3 before tackling QG^. The 
leading candidate for a QG3 theory is the Boulatov model with gauge group SU(2). Its amplitudes arc those of the 
Ponzano-Regge model [4l[ and like matrix models it is topological. The recent identification of a class of graphs 
which for 3d tensorial models generalizes 2d planar graphs, in the sense of having the highest superficial degree of 
divergence (4^ . is an important step for this program. With this paper, we address the question of power counting 
from a slightly different perspective, focusing on general bounds rather than explicit estimates for particular classes 
of graphs. 



D. Constructive field theory 



Constructive field theory was initially the name of a mathematical physics program launched forty years ago by 
Arthur Wightman to define rigorously particular field theory models in increasing order of difficulty, and check that 
they obeyed the "Wightman axioms" [43j . 

The philosophy was to introduce as many cutoffs as necessary for the quantum field correlation functions of these 
models to be well defined, and then to develop the necessary methods to lift these cutoffs. But ordinary perturbation 
theory could not be used directly since for Bosonic field theories it diverges, i.e. it has zero radius of convergence. 
In constructive theory, the divergence of perturbation series is addressed with the same care as the UV divergence of 
individual Feynman amplitudes. Techniques had to be developed to rewrite the formal expansions of ordinary quantum 
field theory as convergent series. Certainly the renormalization group of Wilson is the most powerful such tool. It had 
to be adapted to constructive purposes [44| and was rewritten in the form of rigorous multiscale expansions, slowly 
improved over the years in many ways, and generalized to new situations, such as the RG around the Fermi surface 
of condensed matter [4^ . 

Early constructive field theory succeeded in building rigorously super-renormalizable field theories such as the 
emblematic 02 or models. It also elucidated their relationship to perturbation theory: the Schwinger functions of 
these models are the Borel sum of their perturbation theory [1^, ll^. But 01 itself could not be built, since its coupling 



^ The metric system uses M = 10 for reasons linked to human anatomy, but computers might prefer M = 2 and mathematicians M = e. 

^ Such a transition is usually called geometrogenesis. But even a QG4 ultraviolet topological fixed point might not be the end of physics. It 
could hide at still higher (transplanckian?) phase transition which we could e.g. nickname topologicogenesis, perhaps as a condensation 
of a combinatoric phase. Quantum graphity is a candidate catchy name for such a combinatoric phase 40]. But such a highly speculative 
scenario presumably cannot be explored with QG4 alone. In our notation, graphity or such meta-theories of gravity could be nicknamed 
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constant does not remain small in the ultraviolet regime. Being asymptotically free, non-Abelian gauge theories do 
not have this problem. Nevertheless, although some partial results were obtained, they could not be built in the 
full constructive sense either, due to technical difficulties such as Gribov ambiguities; neither could the interesting 
infrared confining regime of the theory be understood rigorously. Probably the first 4 dimensional field theory that 
will be built completely through constructive methods will be the Grosse-Wulkenhaar model ^38^ , a non-commutative 
field theory which, ironically, should not satisfy the Wightman axioms of the initial constructive program. 

The constructive program is largely unknown in the quantum gravity community, with the notable exception of 
[4^ . However, we think that constructive theory embodies a deeper point of view on QFT than the usual one. A 
modern constructive technique such as the cactus expansion used in this paper allows to resum perturbation theory 
by reorganizing it in a precise, explicit though completely different manner than the usual Feynman graphs. Rather 
than resumming families of entire Feynman amplitudes, this constructive recipe splits Feynman amplitudes into lots 
of finer pieces, according to a precise rule involving new interpolating parameters, and resum infinite families of these 
finer pieces. In the end the cactus expansion itself is a sum of new "amplitudes" indexed by trees or "cacti". The 
series obtained in this way converges absolutely. We believe that this convergent expansion (which is no longer a power 
series in the coupling constant) should be considered the mathematically correct version of perturbation theory.'^ 

Introducing these constructive methods for more exotic theories such as the GFT's should bring further constraints 
on the mathematical consistency of these theories. These additional constraints are precious in view of our lack of 
experiments to test the various theories of quantum gravity. 



E. Plan of the paper 

The plan of our paper is as follows. Section II introduces the Boulatov and Boulatov-Freidel-Louapre (BFL) models. 
In section III we establish a set of uniform perturbative bounds, and exhibit classes of graphs which saturate them, 
showing that the Boulatov and BFL models have very different 'ultraviolet' behavior. In section IV we introduce 
the constructive 'cactus' expansion, which allows to construct the Borel sum of a general 0^ model. In section V, 
we perform this construction for the BFL model, and obtain a bound for the Taylor-Borel remainders of the theory 
which confirms the perturbative estimate. We present our conclusions and future program in section VI. 

Our main results are Theorems IV. II IIII.ll IIII.2I and IIII.31 They establish rigorously that the Boulatov model 
(without mass renormalization) scales perturbatively as h?^"^ per vertex at large cutoff A and that the BFL model 
scales both perturbatively and constructively as A"^ per vertex at large cutoff A - which is much worse. 

The main technique for both perturbative and constructive bounds used in this paper is the Cauchy-Schwarz 
inequality, applied to appropriate cuts of the lines and/or the vertices of the graph. This is also the main tool used in 
[49] in a study of the volume conjecture for classical spin networks, a program which has some obvious overlap with 
ours. 



II. THE REGULARIZED BOULATOV MODEL 



A. Boulatov's original model 

Boulatov's model (26j is a group field theory whose dynamical variable is a real-valued function (j) over the group 

SU(2)'^. Furthermore, the field is required to be invariant under the right diagonal action of SU(2) and cyclic 
permutations c of its arguments, i.e.: 

(t>{9ih,g2h,g3h) = (I){gi,g2,g3) , 0(gc(i), 3c(2), 5c(3)) = 52, 53)- (H.I) 
The "dynamics" of the field is governed by the following non-local action: 

1 f ^ X f ^ 

SB[(f>] '=2 Yi'^dt <^^(5i,5'2,33) + g / Yidg, (/-(ffi, 52, 53)0(53, 54, 55)'/'(ff5, 52, 56)0(36, 34, ffi), (11.2) 



From this perspective, one should not confuse the words "constructive" and "non-perturbative" . In field theory the word non-perturbative 
usually refers to effects which cannot be seen in perturbation theory, such as instanton effects, or to the simulation of functional integrals 
e.g. through Monte-Carlo computations. In contrast, constructive theory is a convergent reorganization of perturbation theory and is 
not yet a theory of non-perturbative effects, which, to this day, remains to be mathematically developed. 
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where dgi denotes the SU(2) Haar measure. Note that the quadratic term is a kind of pure mass term with no "kinetic" 
component. In the quartic term, denoted T[(j3\ below, the six integration variables are repeated twice, following the 
pattern of the edges of a tetrahedron, as illustrated in Fig. [TJ (The factor 8 is only for simpler final formulas.) 

In the language of degenerate Gaussian measures (Appendix A), the Boulatov partition function can be formally 
defined as 

ZsiX) ■■= J d^M e-^^[^]/8, (II.3) 
where dfj,[(j)\ is the normalized Gaussian measure whose covariance is given by the symmetrizer 

,9c(2)h,gc{3)h)- (11.4) 
Thus, in (|!L3|1 . the field (/) is generic, and the invariances (jll.l|) are implemented by the orthogonal projector C^. 



B. A positive regularization 

The partition function pi.3p is only formal and needs regularizations to become mathematically well-defined. 
The problem is twofold: the Fourier space of the field is non-compact, although discrete, and hence 'ultraviolet' 
divergences arise; Boulatov's quartic interaction T[(j}\ is not positive, or unstable, and is thus unsuited for constructive 
considerations. 

To cure the first problem, we follow J48j and introduce an ultra-violet cutoff A truncating the Peter- Weyl (or 
Fourier) decomposition of the field: 

A 

'^(51,52,53)= E tH'^n.j..j.D^'i9i)D'H92)D^H93))- (11.5) 

ii,i2 j3 

In this formula, the sum runs over the spins ji, j2,i3 up to A; D^{g) denotes the (2j -I- l)-dimensional matrix 
representation of g\ ^j^j^j^ are the Fourier modes of the field (j) viewed as complex- valued tensors and tr denotes the 
trace in the space carrying the tensor product representation associated to the spins ji, j2, J3. In the following, we 
denote TiS^'^ the subspace of i^(SU(2)'^) resulting from this truncation, and Ti.'l^^ := Ti*-^-* n Im C^. The number of 
degrees of freedom left is thus given by 



diuiH.i^'^ ^0{K^). (II.6) 

Now the second problem. That Boulatov's interaction T[(j)\ is not positive can be seen from its Fourier space 
formulation where the interaction term reduces to an oscillatory {6j} symbol [26| . To fix this shortcoming, Freidel 
and Louapre propose to add the following 'pillow' term^ to the action [48| 

6 

PW\ — / n^^^ </'(ffi, 52, 53)'/>(53, 54, 55)0(55,34, 36)0(56, 52, 5i)- (II-7) 

i—l 

Indeed, they show that when \5\ < 1, Is[4'] ■— P[4'] + is positive. To this aim, they introduce the 'squaring' 

operator S mapping to the function 5*0 on SU(2)* defined by 5*0(51,52,53,54) := / dg 0(51,52,5)0(5,53,54)- In 
terms of this new field, the modified interaction Is reads 

/5[0] = (501(1 -f5r)50)4, (n.8) 

where T is the involution transposing the central arguments of 50: T50(5i, 52, 53, 54) '■— 50(51,53,52,54), and (•|-)4 
is the standard inner product in I/^(SU(2)^). (1 -|- 5T) being a positive operator, the modified quartic interaction Is 
is clearly positive. 



^ The word 'pillow' refers to the geometric interpretation of the GFT vertex: if Boulatov's T[(j>] is a tetrahedron, then Freidel and Louapre's 
P[<j>] are two tetrahedra glued along two triangles - a pillow. 
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Covariance Tetrahedron Pillow 

FIG. 1: The covariance C^, the Boulatov tetrahedral vertex T[(ji], and the Freidel-Louapre pillow P [(/>]. The labels on the 
vertices match the ordering of the group elements in the integrand of T [</!>] and P[(j>]- 

Combining the cutoff on spins A and the Boulatov-Freidel-Louapre (BFL) interaction, we get the, now well-defined, 
regularized partition function 

4fUa):= (II.9) 

In [i^, Freidel and Louapre showed that Z^-^ is Borel summable. Their method, however, does not allow to prove 
Borel summability of the (more physically relevant) free energy ■— ^^S^bfl^ provides no control over the 

A —> oo limit. In this paper, we construct J'l^j^ explicitly, and determine the scaling behaviour of its Borel radius 
as A — > cx). But before facing these constructive considerations, let us pause a moment to establish the perturbative 
power counting of the Boulatov and BFL models. 

III. PERTURBATIVE BOUNDS 

Let X"Ag denote the amplitude of an order-n Feynman graph G for these models. In this section we prove the 
following 

Theorem III.l There exists a constant K such that for any connected Boulatov vacuum graph G of order n 

l^il < i^"A^+3"/2. (IILl) 

Theorem III. 2 There exists a constant K such that for any connected BFL vacuum graph G of order n 

l^lg^^l < X"AS+3". (IIL2) 

Theorem III. 3 These hounds are optimal in the sense that there exist graphs Gn and G'j of order n such that 

Here, we focus only on vacuum graphs, but these results can easily be extended to the general case. 

Let G be an amputated 0^ graph, i.e. a graph with coordination 4 at each inner vertex and coordination less than 
4 at its outer vertices. We say that a set A of vertices of G is connected if the subgraph made of these vertices and 
all their inner lines (that is all lines of G starting and ending at a vertex of A) is connected. We collect now some 
graph theoretic definitions that will be useful below. 

Definition III.l A generalized tadpole is a two-point graph with only one external vertex. 

Definition III. 2 An (A, B)-cut of a two-point connected graph G with two external vertices va and vb is a partition 
of the vertices of G into two subsets A and B such that va € A, vb € B and A and B are connected. 

Definition III. 3 A line joining a vertex of A to a vertex of B in the graph is called a frontier line for the {A, B)-cut. 
A vertex of B is called a frontier vertex with respect to the cut if there is a frontier line attached to that vertex. 

Remark by parity that there must be an odd number of such frontier lines. 
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Deflnition III. 4 An exhausting sequence of cuts for a connected two-point graph G of order n is a sequence Aq — 
^ Ai ^ A2 ^ • • • ^ An-i ^ An = G such that {Ap, Bp := G \ Ap) is a cut of G for any p = 1, ■ ■ ■ , n — 1. 




Hence an exhausting sequence of cuts is a total ordering of the vertices of G, such that each vertex is roughly speaking 
pulled successively through the frontier without disconnecting A nor B. 

Lemma III.l If G is a connected two-point graph which has no generalized tadpoles subgraphs S C G, there exists 
an exhausting sequence of cuts for G. 

Proof We proceed by induction. We suppose such a sequence Aq = ^ Ai ^ A2 ■ ■ ■ Ap has been built for 
< p < n — 1. We want to find a frontier vertex Wp+i, such that Ap+i — ApU {vp+i}. We consider a rooted tree Tp 
spanning Bp, hence with {n — p) — 1 lines, with root vb (see Fig. 

This rooted tree induces a partial ordering of the vertices of Bp. Since Bp is finite there exists a maximal frontier 
vertex v„iax with respect to that ordering, that is a frontier vertex such that the "branch above Vmax" in Tp does not 
contain any other frontier vertex. 

We notice first that Vmax 7^ vb- Indeed, otherwise it would mean vb is the only frontier vertex left in Bp. The 
number of frontier lines being odd, either vb would have three frontier lines which would mean Bp — {u_b} and would 
contradict p < n — 1, or it would have one frontier line and there would be a generalized tadpole aX vb- 

Second, we clHim that picking ^p-f i — "^-max 

is a valid choice, namely that [Ap U {wp+i}, Bp \ {vp+i}) is a cut. It is 
clear that Ap U {wp+i} is connected. So it remains to check that -Bp+i = Bp \ {up+i} is still connected through its 
inner lines. Call (.p+i the unique line hooked to Wp+i in the path in Tp from Vp+i to the root vb- Cutting splits 
Tp into two connected components, one of which, i?p, contains Vb, and the other is a rooted tree Sp with root and 
only frontier vertex Wp+i. Since Up+i is a frontier vertex, it has at most three lines in Bp, hence there are at most 
two lines from Vp+i to distinct from ^p+i. The tree Rp does not contain any line hooked to tip+i, hence its lines 

remain inner lines of i?p+i so all the vertices in Rp remain in a single connected component of Sp+i. If Bp+i is not 
connected, this would mean, first, that Sp must contain other vertices than Wp+i and, second, that removing in Sp its 
root Vp+i and the (at most two) lines of Sp hooked to it, we would obtain one or two connected components, made 
of the vertices of Sp \ {up+i} plus their inner lines, which no longer hook to Rp through inner lines of i?p+i. Since 
these components have no frontier vertices hence no frontier lines, and since there are no 1-point subgraphs in the 0^ 
theory, it would mean that this component is in fact unique and must have been hooked to G through exactly two 
lines from Wp+i to -Bp+i distinct from ^p+i, hence it would have been a generalized tadpole. □ 



FIG. 2: {Ap, _Bp)-cut in a two-point graph: the dotted lines are the frontier lines, and the solid lines in Bp represent the spanning 
tree Tp, the remaining lines (loop lines) being omitted. When the thinner line lp+\ is deleted, Tp splits into Rp (bottom) and 
Sp (top). 

Now, the 'tetrahedral' and 'pillow' vertices and the covariance can be combined into operators 






(III.3) 



where a = 1, 2, 3 is a channel index, see Fig. [3l in such a way that any BFL Feynman amplitude is obtained as the 
trace of a certain product of these operators and their adjoints. They satisfy the following norm bounds: 
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Lemma III. 2 

||T2,2|| < K, ||ri,3|| < i^A3/2 

\\Pl,\\ < Kk\ \\Pl4 < K, 

11^1,211 ||Pi,3||<ifA3. (IIL4) 





FIG. 3: The vertex operators. In the pillow case, the upper index labels the three different channels. 
Proof To evaluate these norms, we can use the formula 

l/2n 



\H\\ = lim Tt[HH*] 



(rn.5) 



In the Boulatov model, this leads respectively to the computation of the chains G„ for T2 2 and G'^ for Ti 3 as in 

Fig. m 



G„ 



G'n = c QQQQQQQ 



FIG. 4: The chains of graphs Gn and G'„, each with n vertices. 
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These chain are very simple and can be evaluated easily directly; they are also particular cases of 'type I' graphs, 
whose amplitudes were evaluated in [i^. In the Boulatov case, one finds A^^ ~ if", A^, ~ is:«A'^+3«/2 (the factor 
6 comes from the final trace). 

In the BFL case, the dominant graphs are the chains G„ in the channel P22 and with Pi, 3, for which Aq^^'^ ~ 
^ji^6+3ji^ The other chains are subdominant, and yield the results given in pil.2p . □ 

Remarking that a two-point graph is a Ti.^^ —i- Ti-Q^^ operator, this shows that for a two-point graph without 
tadpoles, the norm of the corresponding operator is bounded by X"A'^"/^ and i^"A'^" in the Boulatov and BFL cases 
respectively. Indeed, the existence of an exhausting sequence of cuts precisely allows to write the operator for that 
two-point graph as the composition of n operators, each made of a tensor product of certain number of times the 
identity times a single vertex operator. Applying the bounds of (jIII.2[) . we conclude that Theorems IIII.ll and IIII.2I 
must hold for any graphs without generalized tadpoles: a vacuum graph is the trace of a two-point subgraph, and 
that trace costs at most A^. 

But a two-point generalized tadpole graph G with sub-generalized tadpoles is again a two-point subgraph without 
tadpole joined to a T2.2 or P2,2 operator, so its norm can be bounded again in the same way by an easy induction. 
This completes the proof of Theorems IIII.ll and IIII.2I Finally the examples of graphs G„ and GJ^ prove Theorem 

MM 

The Theorem IIII.2I which express the perturbative scaling behaviour of the BFL model, can be upgraded to the 
constructive level - which is what we now turn to. 



IV. CONSTRUCTIVE FIELD THEORY VIA CACTUS EXPANSIONS 



The historic method in Bosonic constructive field theory is to first introduce a discretization of space-time, e.g. 
through a lattice of cubes, then test the couplings between the corresponding functional integrals restricted to these 
cubes. Since the interaction is local these couplings occur only through propagators of the theory joining different 
cubes. This expansion, called the cluster expansion, results in the theory being written as a polymer gas with hardcore 
constraints. For that gas to be dilute at small coupling, the normalization of the free functional integrals must be 
factored out. Finally the connected functions are computed by expanding away the hardcore constraint through a 
so-called Mayer expansion 0, [s^, [5l|, [13] ■ This last step implies a replica trick in some way or another. Both steps are 
essential; organizing functional integrals around decoupled degrees of freedom is the essence of perturbation theory 
and computing connected functions, or logarithms of the partition function is also a key step in quantum field theory. 
The renormalization group can be seen as essentially the iteration of these two steps, namely functional integration 
over a fluctuation field followed by the computation of a logarithm in order to find an effective action. 

However this historic constructive method, at least in its standard formulation, is unsuitcd for non-local field 
theories. The discretization of the space 811(2)"^ into hypercubes does not make sense, and the non-local interaction 
won't factorize anyway. A similar difficulty arose in the Grosse-Wulkenhaar non-commutative field theory [sl], which 
is non-local in direct space or a matrix model in the matrix base. It was in order to solve this problem that another 
constructive method was recently invented, the cactus (or 'loop vertex') expansion [53]. This method gives correct 
estimates for matrix-like models, and also (as shown in this paper) can be adapted to more general models of the 
tensorial type such as the (regularized) Boulatov model. 

The cactus expansion can be thought as a completely explicit convergent reorganization of standard perturbation 
theory. Its beauty is that, being closer to Feynman graphs than the cluster and Mayer expansions, it is better suited 
to treat more exotic situations such as the non-local interactions of group field theory. Another advantage of the 
cactus expansion is that in some sense it performs the two essential steps of constructive theory (the cluster and 
Mayer expansions) at once: being expressed as a sum of graph amplitudes factorizing over connected components, 
the computation of the connected functions (such as the free energy) is straightforward - one simply restricts the sum 
to connected graphs. 

We illustrate below this cactus expansion on a simple toy model, namely the cj)'^ field theory in dimension.^ 



For a historical perspective on constructive methods applied to this 0-dimensional model, see [s^ . 
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A. Four ingredients for a cactus expansion 

The recipe to perform the cactus expansion of a free function (or any connected correlation function) has four 
ingredients: (1) an intermediate field, (2) a resolvent bound, (3) a replica trick and (4) a forest formula. 

1. The intermediate field representation is a well-known trick to represent a quartic interaction j (p^ in terms of a 
cubic one J (j)a4>: 



e 



where ct is a field with Gaussian ultralocal measure dv (which in zero dimension would simply be the Gaussian 
measure over M). Performing the Gaussian integration over one can then express the original path integral 
Z{\) in terms of a only, for instance in zero dimension as 

Z{\) = / dv{a) I dti{cj)) e-5»VA/0^0 ^ f dv{cj) Aet{l + iVXa)-^'^ - / dv[a)e-^^ Log(i+iv^a)^ (jy 2) 



In more than zero dimension, a propagator C^^^ would typically sandwich the a field on both sides. This process 
generate a 'loop vertex' [5^ for the intermediate field a of the form Tr Log(l -f i^fXa) where Tr stands for the 
operator trace. This trick can be generalized to any correlation functions and to more complicated models. 

2. The replica trick relies on the properties of degenerate Gaussian measures. Let dv{(j) denote the standard 
Gaussian measure on M, and V G L"(M). Now, 'replicate' n times the variable ct, and consider the degenerate 
(normalized, centered) Gaussian measure dv^ on R" with covariance Cij — {piOj) — 1. The replica trick is the 
statement that^ 

/ dv[G) F(a)" = / di.„(ai, . . . , a„) JJ V{a,). (IV.3) 

A general discussion on degenerate Gaussian measures is in Appendix A. 

3. Consider a smooth function H of '^^■'^~^'> variables h = (/i;), living on the lines I of the complete graph over 
n vertices. The so-called Brydges-Kennedy Taylor forest formula jssl . [56j is a Taylor interpolation of H with 
integral remainders indexed by labeled forests over n vertices (see Appendix B): 

"<^'%§.(nH(ni)H<o. 

In this expression, Tn denotes the set of forests over n vertices, the products are over lines I of each forest F, 
and hF is the "^"^~^^ -uple defined by hf :— miup hp, where p runs over the unique path in F connecting the 
source and target vertices of /. (If they are not connected by F, then hf := 0.) 

One can easily check that for n — 2, this is nothing but the fundamental theorem of calculus: H{1) = ^^(0) + 
/q dh H'{h). For higher values of n, on the other hand, the outcome of pV.4|) is genuinely non-trivial, as the 
case n = 3 already demonstrates (Fig. [5]): 

(1,1,1) = i7(0,0,0)+ ( dhi diH{hi,0,0)+ [ dh2 d2H{0,h2,0) + [ 53i?(0, 0, /ig) 
Jq Jq Jo 

+ / dhi / dh2 df2H{hi,h2,m.iii{hi,h2)) + / dhi / dh^ df^H {hi, m.iii{hi, hs), hs) 
Jo Jo Jo Jo 

1 .1 

dh2 dhs disH{mm(h2,h3),h2,h3). (IV.5) 



Jo 



7 



Let us emphasize that this is not a form of Fubini's theorem, which expresses an n-dimensional integral as a product of n integrals. 
Here, we replace one 1-dimensional integral by an n-dimensional one, with {n — 1) delta functions. 
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.-. /. A A A 

FIG. 5: The complete graph over 3 vertices, and its 7 forests, matching the 7 terms in (|IV.5|l . 

4. If A is a complex number with positive real part, and a a real number, one has 

|l + i\/A(7|-i < \/2. (IV.6) 

This inequality easily extends to the case where cr is replaced by a Hermitian matrix S. Indeed, since the 
spectral radius of a normal matrix (and the resolvent (1 + i\fX£j)~^ is normal if E is Hermitian) equals its 
operator norm, taking the supremum over the spectrum of E in (jIV.6|) gives the 'resolvent bound' 



||(l + ^VAI])-l|| < V2, (IV.7) 

in which || • || denotes the operator norm. As we shall see in the following, this bound is key to the summability 
of the cactus expansion. 



B. A toy example 



As a toy example, consider the 0-dimensional 0'* partition function, written in terms of the intermediate 'field' a, 

Z(A) := / e-^^'l^ = / dv{a) e^^("\ (IV.8) 

where Va(o') = — ^Log(l + i\f\o^ and dv(a) is the standard Gaussian measure over the real line. The function 
F(\) := logZ(A) thus defined is analytic in the cut plane C\R_, and hence admits, at best, a Borel expansion about 
A = 0. Using the four ingredients presented above, we now show that this is indeed the case. 

Expanding the exponential in powers and swapping integration and summation* and applying the replica trick to 
the order-rt term yields 

Next, consider the modified covariance , parametrized by an "'■""^^ -uple h, defined by C/j := 1 and := /i;, 
where / = \if\ (i ^ j), and let dv'^{ai, . . . ,an,) be the corresponding (normalized, centered) Gaussian measure. 
Applying the forest formula to the function 



/n 



(IV.IO) 
we obtain 



Of course, it is precisely such an interchange between integration and summation that yields the divergent perturbative series. Note that 
here, however, the process is licit because J dfi{a) e^/2\^°s{^+t^<^)\ < go, and so Lebesguc's dominated convergence theorem applies. 
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^W = E^ E (n fdh) (liw) I d^n{0l,...,^n)i{V^{0.). (IV.U) 

n=o ' Fe;F„ V/e.F'^" / VieF v ti=i 

As announced, the summand factorizes along connected components of each forest, yielding a similar formula for 
the free energy i^(A), except for the fact that only connected forests - trees - contribute: 

^w=E;J[ E fn fA fn^l /rf-r(ai,...,ajnv^AM. (iv.12) 

n=l ' TgT„ ViGT-^O / \leT V v=l 

It is this tree T over loop vertices V\{(jy) which we coin a 'cactus'. 

Since the dependence of the covariance of dv^ (cti, . . . , an) in the h variables is linear, applying the derivative 
in IIV.12l is easy using (jA.2P : it amounts to an additional insertions of ^^^^^ in the integral, where s(/) and t{l) 
are respectively the starting and ending vertices of the line I. Hence 

(n^l) /.-.f K.-.'jni'.Ko =/^.r(. (n ny^M- cv.u) 

Consider the loop vertex V\{av), with coordination ky in the tree T. Thanks to the resolvent bound, we have 



\-^Vxiay)\ = {K - 1)!|A|*|1 + iVAa,r*'-" < 2^(fc„ - 1)!|A|*, (IV.14) 



and thus, since there are n — 1 lines in a tree over n vertices, 



Q2 \ " " 



This bound goes through the normalized integrals over the cr's and the ft,'s, and using Cayley's formula for the number 
of trees over n labeled vertices (jB.2|) . we find that the summand in (jIV.12p is bounded by 2"~^|A|"~^. This shows that 
the cactus expansion (jIV.12p of F converges uniformly in a half-disk Dn = {A € C, 5RA > 0, |A| < R}, with R < ^. 

Let us emphasize here that this result does not contradict the field-theoretic wisdom, which goes back to Dyson 
[13], that a typical QFT perturbative series is bound to diverge: the lemma that convergence at one point implies 
analyticity in a disk around the origin, on which Dyson's argument relies, only applies to power series - and the cactus 
expansion (ITO]) is not a power series. 

Of course, one might wonder what the convergence of the cactus expansion teaches us about the perturbative 
expansion of F{X) 



F(A)^^(-A)% (IV.16) 

p>0 

which, in fine, is the one dealt with in standard QFT. The answer is simple: the cactus expansion proves Borel 
summability of the perturbative series, and gives an explicit expression for its Borel sum (Appendix C). Indeed, since 
V\ is analytic in the cut plane, and F is the sum of a power series in V\ converging uniformly in the half-disk Dn, 
we now know that F is analytic in a Nevanlinna-Sokal disk A'loreover, the order-r Taylor-Lagrange remainder 

TrF{X) :— F{X) — X]p=o ^-^'^'(0) can easily be shown to satisfy the Nevanlinna-Sokal criterion (|C.2|) . Indeed, 
consider a cactus amplitude Ct with n loop vertices: 

■■= (U I' ^'^'^ (n ^) / ^^rf ('^i, ■ • ■ , ^n) n Vxia,). (IV.17) 

By (jIV.13[) - pV.14[) . each such amplitude is made of an explicit factor A"^^ times an integral over dv'^'^ {oi, . . . , (t„) 
of a product of 2r7, — 2 resolvents i?;^, (A, cr„) := (1 + i\/ACT„)^^, where Iv denote a half-line hooked to a vertex v. Hence 
for r < 71 — 1, TtCt = Ct, and for r > n 

TXt ^ X'^-^Tr-n+i m / dhA f di^f{au...,an) n i?;„(A,a„). (IV.18) 
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But by (|A.2[) and since F — nf"=i^ Ri^ is solely a function of \f\a we have 

. 2n-2 ^ 2n-2 

diyfiai,...,an) n i?i„(A,fT,) = 6^^^^" ^ [] i?;„(l,'T) U„=o- (IV.19) 

( 1 ( 1 



Hence the Taylor-Lagrange formula applies to the exponential: 



Tke^H = Cdt ii-i^A^i/^e^*^ (IV.20) 



(fc-1) 



and the operator iJ'^ creates exactly k additional insertions of lines t; ^ . The combinatorics of 2fc derivations 

f2i — 1)V 

on a product of 2n — 2 resolvents costs a factor (2rt-i)! times a new Gaussian integral with 2rt — 2 + 2fc resolvents of 
the same type Ri . The A factor can be then transferred back to a \/A factor in the resolvents: 



2n-2+2fc 2n-2+2fc 
1 \+ a r^h'^ a -w—r , ^ , 1 + S r'^^ a -r-r , ^ , , , 

e2>*^c ^ -Q i^^^^^g^t^c ^ -Q i?;,(A,a„) 1.^=0. (IV.21) 

Hence since k = r — n + 1 the convergent series X^^i ;7T X^tsT ^'-^t" '^^^ be bounded exactly as before, except for 
two facts: each term contains a factor A'' and we have also to add to the bounds a factor dt ^ ^f^*^^ ^ | ^(2n'2iy, ~ 
7 '•^''iTi/y' — TTT- This last factor is maximal for the trivial tree with n — \, and certainly bounded by 2*" r!, from which 

(r— n+l)!(2n— 1)! ' ^ j 5 

it follows that for some constant K 

\TrF{X)\ < ifV!|Ar. (IV.22) 

To summarize, the cactus expansion allows not only to trade the asymptotic perturbative series PV.16P for the 
convergent expression pV.12|) . but also to check the Sokal-Nevanlinna criteria, proving Borel summability of the 
former. Of course, in this toy example, Borel summability of Z is obvious and Borel summability of -F = log Z could 
be shown by more elementary methods; the power of the cactus expansion becomes manifest when it comes to the 
constructive analysis of the 01 field theory [s^ , and of the matrix 0* model [53^ . Our aim in this paper is to extend 
its scope to tensor models such as Boulatov's GFT. 

V. CONSTRUCTION OF THE BFL MODEL 

A. Intermediate field representation 

Let us now construct the cactus expansion of the BFL model. Following the recipe explained above, we first 
introduce a ultralocal intermediate field a on SU(2)^ slicing the BFL (p^ vertices into two (jPa vertices as in Fig. [SI 



g-A/40]/8 



j dyf\a]e-i^'^^*\-^-. (V.l) 

Note that, in this intermediate field picture, the tetrahedral and pillow interactions are encapsulated in the ultralocal 
Gaussian measure dvg'^^ through its covariance Ca :— (1 + 5T). 
Introducing the operator E coupling to cr 

^0(51,52,53) := y" dg4:dg5 a{gi,g2,g4,g5)(t){g3,94,95) (V.2) 

in such a way that 

(50|a)4 = (0|E0)3, (V.3) 
where (■|-)3 is the standard inner product in L^(SU(2)^), we obtain after integration over the original field 



-7(A) 

^BFL 



(A)- J e^^M, (V.4) 
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Tetrahedron 



FIG. 6: Slicing the BFL vertices with an intermediate field cr over SU(2)*: the dashed hnes are combined in the covariance Ca ■ 
where the 'loop vertex' is given by 

Vx[<y] := -^Tr Log(l + zVAC^SC^). (V.5) 

One easily checks that S :— C^Y,C^ is a Hermitian operator and therefore that the resolvent bound applies to the 
derivatives of this 'loop vertex' just like in the toy example. 



B. Cactus expansion 



Following the same steps as in sec. IV yields the cactus expansion of the BFL free energy: 



n=i ■ Ter„ \;gt"^" / 

At this stage, the only difference with the 0-dimensional case is the insertion of a covariance Ca-{gs{i)] 9t{i)) on each 
line I of the tree, and the integration with respect to the corresponding 4-uples of group elements 5'" := 
attached to the half- lines ly. 

Computing the effect of the fc„ derivatives on the loop vertex VA[f7i,], labeled by the half- lines 1^ connecting it to 
the tree, we obtain 

(n - (n a + . (v^^) 

and the cactus amplitude is given by the following product of traces connected by ultralocal covariances: 



Ct(A) := \\\ dhi] du'^ 



n / dV^'^dV'^ n Tr (1 -HzVAE(a.))-i^^^^ . (V.8) 



C. Cauchy-Schwarz inequalities 

To get some insight into this cactus amplitude, it is handy to introduce a 'dual' representation of a tree T, as 
a planar partition of the disk. The boundary of the disk is obtained by turning around T, while the dotted lines 
partitioning it cross the boundary twice and each line of T exactly once, without crossing each other, see Fig. [7l 

In such a picture, the resolvents (1 -I- iV\Tj)~^ are attached to the arcs on the boundary of the disk, while the 
covariances Ca are attached to the dotted lines. To bound (|V.8[) . we can apply the Cauchy-Schwarz inequality along 
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FIG. 7: The planar representation of a tree. 

a line splitting the disk in two parts with the same number of consecutive resolvents. Indeed, the number of half-lines 
of a tree being even, it is always possible to pick two arcs with resolvents i?i and i?2, and express ()V.8|) as the inner 
product {A\Ri (g) R2\B), where A and B contain the same number of arcs, and thus of resolvents (see Fig. [8]). By the 
Cauchy-Schwarz inequality, we have 



\{A\R, ^ R2\B)\ < ||i?i|||i?2|iyPIiyyMBy<2v/pIiyv/MBy- (V.9) 



A A B* 




B A* 



FIG. 8; Splitting the disk in two parts to apply the Cauchy-Schwarz inequality. On the LHS, black squares are resolvents, 
dashed lines are covariances Co-, and the thick line expresses the amplitude (|V.8P as the inner product between the upper and 
lower parts A and B. On the right, white squares are Hermitian conjugates of resolvents and dotted lines are covariances C^. 

This process has two effects: thanks to the resolvent bound, it trades the original tree for two trees, with the same 
number of vertices, but with two resolvents replaced by the identity (||i?i|| ||i?2|l < 2); each covariance Co- which is 
sandwiched in the inner product is replaced by a modified covariance := 1 + 5T' , where T' identifies the two 
central arguments of <t instead of twisting them (Fig. 

T'{gi, ...,gA;g[,.. . ,54) 5{gigi'^)5{9Ag^'^)5{g2g^^)5{g'2g'i^)- (V.IO) 

We can iterate the process [n — 1) times, until all resolvents are removed, and all covariances replaced by C'„. 
We are then left with a perturbative BFL graph, whose vertices are all -P2^2- 

|Ct(A)| < 2"-i sup^,g^„ (V.U) 

Now, from Theorem |IIL2] we have snpj^^r^ \^t^^\ < K"'A^+^'\ hence |Ct(A)| < 2"-iii:"Af^+3". 

Before we can conclude, we should check that this estimate also holds for the trivial tree with just one vertex {n = 1): 
since it contains no half-line, the loop vertex Log(l -I- za/AE) is not acted upon by a u-derivative, and therefore the 
resolvent bound does not apply. However, using standard convexity inequalities, we have |Log(l-|-z)| < i-f i(|zp-|-47r^) 



16 




FIG. 9: Covariances Co- (dashed lines) sandwiched in the inner product in the Cauchy-Schwarz inequality (thick line) are 
replaced by modified covariances (dotted lines). 



for z e C, hence 

|lLog(l + ^^/AE)|l <\ + + A^^) <\ + \{Ml + 4^^), (V.12) 

and thus 

|Tr j dvf\a) Log(l + iVAE)| < K AimUi^'^ = 0{K^). (V.13) 

Absorbing the uniform factor dimTYQ^'' — 0{K^) appearing in all estimates in the definition of a free energy per 
degree of freedom ^/bfl — ^^Tat-^bfl> '^^ have proved^ 

Theorem V.l The cactus expansion of the BFL free energy per degree of freedom Ggpi^ is uniformly convergent in 
a half-disk {A G C,5RA > 0, |A| < KA~^}, where it defines the Borel sum of the BFL perturbative series. 



VI. CONCLUSION 



We have proved that individual graphs at order n of perturbation theory for the Boulatov and BFL models are 
uniformly bounded by X'^A^^^'^^^ and /C"A^+^" respectively, and that the Taylor-Borel remainders for the BFL 
model are bounded by n!i^r"A^+^". These bounds are optimal, being saturated by a class of graphs. Since the 
perturbative power counting of the Boulatov model is much better than the one of the BFL model, we conclude that 
the Freidel-Louapre constructive regularization with the 'pillow' term of the Boulatov model is not optimal from the 
perspective of power counting. 

We conjecture that 

Conjecture Any Boulatov graph at order n without two-point subgraph (hence with A'^ > 4 external legs) is bounded 
at fixed external legs arguments by C"A"^^. 

Remark that this conjecture is slightly different from the conjectures of |42i|, as the latter paper is more oriented 
towards identifying the leading graphs (called type I) and finding their exact power counting. 

Our future program consists in 

• Identifying a better constructive regularization than the BFL "pillow term" , i.e. such that the Taylor-Borel 
remainders obey the same scaling in ri!C"A^+'^"/^. 

• Proving the conjecture, presumably through a multiscale analysis, and identify a still better constructive reg- 
ularization, i.e such that the Taylor-Borel remainders of e.g. the 4-point function, after two point function 
renormalization obey the same scaling in n!C"A"'~'^. 

• Generalizing in the appropriate way these conjectures, first to the Ooguri model in dimension 4, then to the 
EPR-FK models and investigate their scaling properties in detail through the appropriate multiscale analysis. 



The Nevanlinna-Sokal criteria can be checked exactly as in the toy example. 
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APPENDIX A: DEGENERATE GAUSSIAN MEASURES 

If A is a n X 71 positive definite matrix, a (centered, normahzed) Gaussian measure fi^ over R" can be defined in 
terms of its Radon-Nikodym derivative with respect to the Lebesgue measure A 



-7— (x) = —TTre 2 . (A.l) 

From this definition, it follows that the covariance matrix C — A^^ is positive definite. Wick's theorem then implies 
that all higher moments are completely determined by C. In fact, the measure itself is uniquely characterized by the 
covariance, as an example of a positive solution to the Hamburger moment problem. Since the last two propositions 
(Wick's theorem and reconstruction from the moments) would also hold true if C were positive semi-definite, a more 
general definition of a Gaussian measure is as follows: given a positive semi-definite matrix C, the Gaussian measure 
fic with covariance C is the unique measure whose moments are given by Wick contractions of C. Integration with 
respect to dfic is then really the same as applying the exponential of a quadratic differential operator to the integrand. 
More precisely, for any smooth and summable / 

dA^c(x) /(x)=e^^^^/(x)|_„ (A.2) 

When C has zero eigenvalues, dfic is a 'degenerate Gaussian measure', with a Dirac delta component along each 
zero mode. In such cases, /ic is not absolutely continuous with respect to the Lebesgue measure, and (lA.ip does not 
make sense, while (|A.2|) does. 

A first useful example of such a situation is the replica trick considered in section [TVl which is the key to the rigorous 
construction of connected functions in QFT. It relies on the existence of the degenerate Gaussian measure defined by 
the rank-1 matrix = 1. Its eigenvalues are 0, with multiplicity n — 1, and n, with multiplicity 1 and eigenvector 
(1, . . . , 1). It follows that for any smooth and summable function / in n variables: 

dficix) fix) f da e--"/'f{a, ...,<j). (A.3) 

V 27r Ja 



The interest of such a seemingly trivial measure is the possibility to perturb off diagonal coefficients of this Gaussian 
measure to break the replica symmetry. 

Another instructive example is precisely given by group field theory. It is often formulated as a functional integral 
over invariant fields, that is loosely speaking as 

"D„,,cl^e-^'/^-'^'^\ (A.4) 

where the tensorial nature of the (j>'^~^^ term reflects the gluing of d-simplices along common (d — l)-faces, and the 
Gaussian part looks like a trivial "local" mass term. The nontrivial content of the theory is then completely hidden 
in the "Lebesgue measure" Di„„ ((/>), which restricts functional integration to fields invariant under e.g. diagonal right 
multiplication: 

'/'(ffi, 9d) ^ (j){9ih,...gdh) Vh . (A.5) 

But this loose writing is confusing. Remember that in the case of an ordinary field theory, such as the ordinary 0"* 
theory, the loose writing of functional integration as 

^^g~/i0(-AW)0-A/0^ (A.6) 



A. Connes always remarks, especially when discussing Galois's work, that it is not symmetry which is important, but how to break it. 



18 



is already misleading because there is no well-defined "infinite dimensional Lebesgue measure" Dcj) = d(f>{x) in 
any reasonable sense. As well known in constructive theory, f D(f> and /^(-A+m )</> gj^ould be combined into the 
Gaussian measure d^c with propagator C = (—A + rn?)~^ which is mathematically well-defined as a measure on a 
space of distributions through the Minlos Theorem. One should therefore write 

j dp.c{<l^)e-^i^\ (A.7) 

Then one identifies the problems: the interaction term in dimension two or more is not Radon-Nikodym with 
respect to that Gaussian measure because of ultraviolet divergences, and moreover a thermodynamic limit has also to 
be performed because the infinite volume integral 0* diverges with probability one. The expression (|A.7[) at least 
opens a sensible avenue, which is to apply an ultraviolet cutoff on C in d/ic, start with a finite volume instead of M'*, 
and perform the two limits in the right way (that is step by step, according to the renormalization group). 

In (|A.4[) . Dinv4> is the product of the previously ill-defined "infinite dimensional Lebesgue measure" D<j) by S- 
functions constraints implementing the constraints (jA.5|) . So in a way it has no Radon-Nikodym density, but with 
respect to a measure which is not well defined! This is why the writing (|A.4p is doubly confusing. The way out is 
again to remark that because (S-functions are Gaussian measure, one can combine the constraints, the 'mass' term 
and the Lebesgue measure into a Gaussian measure dfic with propagator 

C{gu.-,gd;g[,-,g'd)^ j dh 5{g^^hg[) ■ ■ ■ 5{g^^hg'^) (A.8) 

which is now well-defined. 



APPENDIX B: TREES AND FORESTS 



For convenience, we collect here some basic definitions of graph theory. We use the quantum field theorists' 
vocabulary, recalling the mathematicians' in parentheses. 

A graph (pseudograph) G is a set of vertices V and of lines (edges) E, together with an incidence relation between 
them, possibly with several lines connecting the same vertices (multiple edges) and lines connecting a vertex to itself 
- tadpoles (loops). A subgraph G' of G is a subset of edges of G, together with the attached vertices. A loop (cycle) 
is a connected subset of n lines and n vertices which cannot be disconnected by removing any line. 

A forest is a graph without loops. A spanning forest of G is a sub-forest of G that contains all the vertices of G. A 
tree is a connected forest, or equivalently a graph with \E\ = \V\ — 1. A vertex with coordination 1 in a tree is a leaf. 
A rooted tree is a tree with a distinguished vertex, its root. 

A useful result in graph theory is Cayley's formula, giving the number r„ of different labeled trees on n vertices, 

Tn = n^-^ (B.l) 
and the number of such trees with fixed coordinations ky 

^ n:..(^.-i)r 

Four proofs of this formula are given in (soj . 



APPENDIX C: BOREL RESUMMATION 



It is a classic theorem of Borel that there is an infinite number of smooth real functions asymptotic to any power 
series. But analytic functions are rigid: within their domain of analyticity, all the information about the function is 
encapsulated in the countable list of its Taylor coefhcients. Ordinary summation provides a one-to-one correspondence 
(at least inside a convergence disk) between a convergent power series and a unique "preferred" function asymptotic 
to that series, namely the analytic one. 

Borel summability is a natural way to extend this one-to-one correspondence between convergent series and analytic 
functions. When an analytic function admits a Taylor series around a point on the boundary of its domain of 
analyticity, this series has zero radius of convergence. But under some conditions, all the information about the 
function can be still be encapsulated by the series. Borel summation picks up the unique "preferred" function 
asymptotic to that series, namely its Borel sum. 
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Borel summability is ubiquitous in theoretical physics, where most expansions have zero radius of con verg ence. 
That it was the case of the perturbative series of quantum field theory was realized by Dyson back in 1952 ,57]. 
The most natural criteria for Borel summability were formulated by Nevanlinna [60'|, and rediscovered by Sokal 

m- 

Theorem C.l (Nevanlinna-Sokal, direct) Let f be analytic in the disk Cr :— {y\Key^^ > 1/-R}. Suppose f 
admits an asymptotic power series ^Qky'^ (its Taylor series at the origin) 

f{y) = Y,aky'' + Rr{y) (C.l) 

fc=0 

such that the hound 

\Rr{y)\<caW\\y\^ (C.2) 

holds uniformly in r and y G Cr, for some constants a and C . Then f is Borel summable, i.e. the power series 

■= Sfc*^fclT converges for \t\ < ^, and admits an analytic continuation in the strip So- '■= {t\ dist (ijIR"*") < 
satisfying the bound 

\B{t)\ < const. eTT for t e M+ (C.3) 
Moreover, f is represented in Cr by the absolutely convergent integral 

f{y) = - / dt e'^B{t) (C.4) 

y Jo 

B is then called the Borel transform of /, and the complex t plane is called the Borel plane. 
There is a reciprocal to this theorem: 

Theorem C.2 (Nevanlinna-Sokal, reciprocal) Consider the power series J^'^kV^- If the power series J^'^kjr 
converges in a disk \t\ < ^, admits an analytic continuation B{t) in the strip S„ and satisfies the bound \C.S^) in this 
strip, then the function f defined by the integral representation |C.^[ j is analytic in Cr, has 'Y^aty^ as Taylor series 
at the origin and satisfies the uniform remainder estimates. 

In this case we say that the series J^'^kv'^ is Borel summable, and call the series X^'^^It i^s Borel transform and 
that the function / its Borel sum. 

In conclusion, Borel summable series and Borel summable functions are in correspondence just like are ordinary series 
and germs of analytic functions. However, the analytic continuation in the Borel strip involved in the construction of 
the function from its series is usually intractable, and so in practice only the direct theorem is used. 

For interactions of higher degree than 0"*, which will be required for group field theory in dimension 4 and higher, 
extended notions of Borel resummation, such as Borel-LeRoy resummation of order a, must be applied. In essence, 
they combine suitable conformal transformations of the standard disk Cr with the standard Nevanlinna-Sokal criteria. 
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